Dynamical band nipping in fermionic lattice systems: An ac-field-driven change of the 

interaction from repulsive to attractive 
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We show theoretically that the sudden application of an appropriate ac field to correlated lattice 
fermions flips the band structure and effectively switches the interaction from repulsive to attractive. 
The nonadiabatically driven system is characterized by a negative temperature with a population 
inversion. We numerically demonstrate the converted interaction in an ac-driven Hubbard model 
with the nonequilibrium dynamical mean-field theory solved by the continuous-time quantum Monte 
Carlo method. Based on this, we propose an efficient ramp-up protocol for ac fields that can 
suppress heating, which leads to an effectively attractive Hubbard model with a temperature below 
the superconducting transition temperature of the equilibrium system. 
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Introduction. — There is an increasing fascination with 
the dynamics of fermions driven by external fields as a re- 
sult of recent developments in time-resolved experimen- 
tal techniques in, e.g., ultracold atom physics [l| and 
electrons in a crystal [2j. There, the external fields are 
employed not only for probing the response against per- 
turbation, but also for creating excited states to control 
phases of the system. In particular, continuously driven 
systems isolated from the environment exhibit nonequi- 
librium statistical distributions that dramatically alter 
their physical properties. One basic, long-known example 
is the negative temperature (T) distribution |3|, where 
higher energy levels are occupied (population inversion) . 
If such distributions are realized in correlated fermion 
systems, it should have a huge impact on many-body 
physics since, as we shall show, the originally repulsive 
interfermion interaction effectively turns into an attrac- 
tive one as a result of the negative T. 

To start with, a realization of negative T requires as 
a crucial condition that the energy spectrum is upper 
bounded 0). A spin system, for which the concept of 
negative T was originally introduced satisfies this 
condition. Another promising candidate for negative T 
is a many-particle system on a lattice, where particles 
form a band structure with a finite bandwidth, as far as 
the interband excitations induced by external fields are 
forbidden or negligible. For instance, Ref. [B| proposes 
that a harmonic potential trapping bosonic atoms in an 
optical lattice can be inverted to realize a negative T. 

In the present Letter, we take a totally different ap- 
proach, namely, a sudden application (quench) of a sinu- 
soidal force (ac field) to a lattice fermion system. The 
effect of ac fields has been discussed previously in the 
context of adiabatic processes (e.g., [6]), where the hop- 
ping amplitude is effectively renormalized Q- Here we 
show that the "nonadiabatic" switch-on of the ac field can 
dynamically invert the band structure (dynamical band 
flipping), generating a negative-T distribution. Based on 



this, we suggest that the ac quench is related to an inter- 
action quench ("ac/[/-quench correspondence"), where 
the interaction can be effectively converted between re- 
pulsive and attractive. This is numerically demonstrated 
for the ac-driven Hubbard model, whose time evolution is 
obtained with the nonequilibrium dynamical mean-field 
theory (DMFT) @, i] solved with the continuous-time 
quantum Monte Carlo (QMC) method fioj j . 

One immediate question may be whether the sudden 
switch-on of the ac field causes a violent heating of the 
system. However, we shall show that the heating as- 
sociated with the sign change of the interaction can be 
suppressed if the ac field is turned on in a suitable way. 
Therefore, if the ac-driven system thcrmalizes, it will 
correspond to an effectively attractive Hubbard model 
whose temperature is low enough that the system can 
possibly accommodate superconductivity. 

Formulation. — We consider a model Hamiltonian, 

H(t) = -JH K + UHi + cos(fti) K ( t ) ■ R i n i° W 
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Here c-„ creates a fermion at 



£j (riff- |) [n 3i -\). „ ja 

site j with spin a for electron systems or pseudospin a 



for cold atoms, Rj is a position vector, jij a 
J(> 0) the hopping energy, U{> 0) the on-site inter- 
action. The third term in Eq. ([1]) represents the ac 
field with amplitude K and frequency O. We take, 
for the DMFT, a d-dimensional hypercubic lattice with 
K(t) = K(t)(l, 1, . . . , 1), and consider a half-filled band. 
We assume that the system is initially in equilibrium at 
temperature T, and the ac field is suddenly switched on 
at t = 0, i.e., K(t) = K6{t) with 9{t) the step function. 

It has been shown [7[, for noninteracting systems, that 
the effect of the ac field is simply a renormalization of 
the hopping energy, 



J -> J eS = MK/n)j, 



(2) 



2 



with J n (z) the nth order Bessel function. A naive expla- 
nation of Eq. ([2]) is that in the presence of the ac field 
the original band dispersion tf. = — 2JX),=i cosfc^ is re- 
placed with a time-averaged Tk — ^ f^^ 1 d^k~A(t) — 
Jo(K/il)ek, where we take a gauge in which the ac field is 
represented by a vector potential A(t) = — K sin(f2i)/f2. 
The renormalization of J can more rigorously be derived 
from the Floquet theory for ac fields [ll|. Curiously, 
Jeff oc Jo(K/fl) then changes sign as A = K/Vl is in- 
creased. This implies that J e g vanishes at Ai = 2.404 . . . 
(see the inset of Fig. |3|) , which is known as dynamical lo- 
calization [7]. Experimentally, the scaling ([2]) was beau- 
tifully observed in a Bose-Einstein condensate (BEC) in 
an optical lattice [l2j ]. 

What happens when Jo{A) < 0? In equilibrium, the 
inverted sign of J does not change the physics, since it 
can be canceled with a gauge transformation cj CT — > — cj CT 
for the j € B sublattice. Out of equilibrium, however, we 
have to consider the distribution of occupied states, and 
the physics can indeed dramatically change as we cross 
the zeros of Jo(A). 

To characterize nonequilibrium distributions, we can 
define a time-resolved spectral function, A{ujJ) = 
-ilm / dt'e iuJt ' J2k G ka{ 1 + T> t ~ t) and an occupation 
N(u,t) = / dt'e^' E fe G< CT (t+|,t-^), in terms of 
the retarded G% a {t,t') = -i6(t-t'){{c ha (t), 4 CT («')}) and 
lesser Green function G^ a (t,t') — i(Cf.(t)ch„(t')). For 
noninteracting fermions we can evaluate [111 ], for 51 3> J, 
the longtime behavior as 



A(u,t) ^ 



\MA)\ Po \\MA)\J' 



N(w,t) A(u,oo)fi 



MA) 



(4) 



where po(u)) — — e fc) ^ s the noninteracting den- 

sity of states, /t(w) = l/(e"/ T + 1) the Fermi distri- 
bution with ks = 1, and we have used the symmetry 
Po(— lj) = po(u>). We notice that the sign of Jo(A) is ir- 
relevant to the spectrum, but it is indeed relevant to the 
occupation. In fact, for Jo (A) < the effective temper- 
ature is negative, since fr(u / Jo(A)) — /_t(w/| Jo(A)\). 
One intuitive way to interpret this is that, as the band 
is flipped when J changes sign due to the ac field, each 
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FIG. 1: A schematic band flipping accompanied with a pop- 
ulation inversion (T < 0) induced by the ac field when 
Jo (A) < and without a population inversion (T > 0). 
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t-tt= \ J (A)\t 
rescaling of time 



//-quench (5) with p\(f) 



| energy reversal 



//-quench (5) with p^{t) 



FIG. 2: The ac/[/-quench correspondence. 



fermion follows the change with no shift in the momen- 
tum (Fig. Q] for T < 0). This is the essential mechanism 
for the negative T distribution that we propose here. 

Interacting systems driven by ac fields have also been 
theoretically studied for the Bose- Hubbard model 0, [HI , 
and the Falicov-Kimball model [III , [lij ]. Let us as- 
sume that J is suddenly quenched by the ac field 
as in Eq. ([2]). Then the density matrix time-evolves 
as pit) = e-M-JosfiK+un^pfn^iti-J^K+uni) _ 
^-mK+u*Hx) p (Q)e±a<.-JHK+u«,Hi) = p ±(;) ; where 

the upper (lower) sign corresponds to Jo(A) > (<)0, and 
time is rescaled as t = \Jo(A)\t. After that, the J quench 
is translated to the interaction quench [l5| as 



U ->■ U eS = U/J (A). 



(5) 



For Jo (A) < 0, the corresponding U quench evolves in 
time with p^(t) where the phase rotates in the reverse 
direction. We can relate p^(t) with the normal time 
evolution p\(t) via an energy reversal; i.e., the signs 
of all the quantities (except for time) that have dimen- 
sion of energy are reversed. For example, energy reversal 
(3) maps kinetic energy E Mn (t) = -i J2 k ^k-A{t)G^ a {t, t) to 
— Eki n (t). We summarize the relation between ac quench 
and U quench in Fig. [2J 

The ac/C/-quench correspondence, which we numeri- 
cally demonstrate in the following, has the intriguing 
consequence that for Jq(A) < the ac field effectively 
switches the many-body interaction from repulsive to at- 
tractive (U c s < 0) ((5]) . In addition, it implies that the 
system thermalizes to a negative-temperature canonical 
distribution: according to the thermalization hypothe- 
the (nonintegrable) Hubbard model evolving 
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with normal p\{t) thermalizes after a U quench in the 
long-time limit with a positive temperature T e ff(> 0). 
Since p~^{t) is related to p\{t) via the energy reversal 
(Fig. [2]), the ac/ [/-quench correspondence implies that 
the system driven by the ac field should finally thermal- 
ize with a negative temperature, — T e g/\ Jo(A)\. 

Numerical result. — To verify the interaction switch- 
ing implied by the ac/C/-quench correspondence in the 
dynamical band-flipping in interacting systems, we have 
numerically computed the time evolution of the system 
using the nonequilibrium DMFT which can treat 
the dynamics of infinite systems with the local correla- 
tions from the many-body interactions fully taken into 
account. In the DMFT calculation, we restrict ourselves 
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FIG. 3: Time evolution of the double occupancy D for various 
values of A = K/Q. (symbols with error bars) for Q. = 2ir, 
U — 1, with the horizontal line indicating the free-electron 
value of D — 0.25. Solid curves represent the corresponding 
result for the U quench with p\ (t) plotted on a rescaled time 
axis t = t/\Jo(A)\. Inset depicts the Bessel function, Jo(A). 



to paramagnetic phases without charge order. The hop- 
ing is scaled with dimension d as J = J*/2Vd (d — > oo) 
. In the following we take J* as the unit of energy and 
write J* as J. Even when fermions are initially weakly 
interacting, the effective interaction ([5]) plunges into the 
strong-coupling regime around the zeros of Jo^A), so that 
we need an impurity solver that is valid for both weak 
and strong interactions. Here we employ the continuous- 
time QMC method based on a weak-coupling expansion 
10], which is numerically exact within statistical errors. 

A simple measure of whether the interaction is re- 
pulsive or attractive is the double occupancy, D{t) = 
(ntf(t)rij±(t)) . Figure[3]displays the time evolution of the 
double occupancy for various amplitudes A = K/fi with 
fixed fi = 2tt, U = 1, and T = 0.1. D(t) starts from the 
equilibrium value, which is, for the repulsive interaction, 
smaller than the noninteracting value, (fijt)( n .j4.) = 0.25. 
For A < A\ (left panel of Fig. [3]), D{t) decreases with 
t, which is natural since the effective interaction ([5]) is 
enhanced. For large effective interaction, the double oc- 
cupancy shows 27r/[/-periodic collapse-and-revival oscil- 
lations [lj| , while the fast oscillations with frequency 2S1 
come from a nonlinear effect of the ac field. 

The double occupancy starts to behave in a dramati- 
cally different manner as we plunge into the Jo (A) < 
regime, i.e. A\ < A < A2 = 5.520 . . . (right panel 
of Fig. [3]): D(t) steeply increases after the ac field is 
switched on, and even goes "beyond" the noninteracting 
value 0.25. This implies that fermions prefer a large dou- 
ble occupancy, evidence that the many-body interaction 




FIG. 4: (a) The momentum distribution, (b) the spectral 
function (thick curve), the occupation (shaded region), and 
(c) the kinetic energy for an ac quench with A — 4, fl = 2n 
and U = 1 (oscillatory curve, along with a time-average in 
black) , as compared with the corresponding U quench (dashed 
curve, along with a sign-inverted plot). The top left inset 
schematically shows a negative-T situation. 



indeed turns into an attraction. When Jq(A) returns to 
positive for A > A2, D(t) again becomes smaller than 
0.25, and the effective interaction goes back to repulsive. 

To endorse the ac/C/-quench correspondence (Fig. [5]) 
quantitatively, we plot D(t) [solid curves in Fig. [31 plot- 
ted as a function of the rescaled time t — t/\Jo(A)\] 
and -Ekin(i) [dashed curve in Fig. EJc)] for the U quench 
calculated with p\{i). For D(t), we can see that the 
[/-quench results agree with those for the (nonoscilla- 
tory components of) the ac quenches to a surprisingly 
good accuracy. The accuracy is especially excellent for 
A > 4, while for A < 3 the small differences seen be- 
tween the U and ac quench come from the fact that 
it takes a few cycles for the ac field to renormalize J. 
For Jo(A) < we again confirm the ac/?7-quench corre- 
spondence in the kinetic energy if we consider the energy 
reversal Ek- m (t) — > —Ek in (t) in the U quench [a sign- 
inverted plot of the dashed curve in Fig. Q|c)] and time- 
average Eki n (t) over one cycle in the ac quench (black 
curve) . 

As for the population inversion, we can also obtain its 
direct evidence by calculating the (gauge-invariant) mo- 
mentum distribution function [ll[ f{k,t) = —iG^(t,t) 

(Mt) = k + A(t)) in Fig. gja), and A(uj,t) and N(oj,t) 
|17j in Fig. BJb) for the same U, f2, and T as in Fig. 02 
with A = 4 [for which J (A) < 0]. We plot /(fe,t) along 
a typical slice [18| in the Brillouin zone. In the initial 
state, the particles distribute around — ^ < k < | with 
a blurred Fermi surface. As the ac field is turned on, 
the distribution is gradually smeared out, but the main 



population continues to stick to — § < k < 



This 



sharply contrasts with what has been observed in BECs 
[12j : there, the system adiabatically follows the lowest 
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FIG. 5: An example of a multistep quench (left panel), and 
the corresponding evolution of T e g(t) (right). The inset indi- 
cates the path for the ac quench on Jo{A). 



energy states with no population inversion (Fig. [T] for 
T > 0) with the peak of f(k, t) shifted by tt in the nipped 
band [19j. However, our results indicate that, with suf- 
ficiently fast ramp up of the ac field in fermion systems, 
the system does not relax to lower energy states, but the 
population can be inverted. The population inversion is 
more directly seen in N(oj,t) [Fig. 0Jb)], where fermions 
tend to occupy higher energy states (oj > 0), justifying 
our picture (Fig.Q]for T < 0). This kind of drastic change 
of the dynamics as a function of the parameter is rem- 
iniscent of phenomena observed in nonlinear mechanics. 
It may thus be interesting to try to connect bifurcation 
theory [iol ] to the band-flipping phenomena found here. 

Suppression of heating — Now that the ac field is 
found to effectively convert a repulsive interaction into 
an attraction, one may contemplate the possibility of ac- 
induced superconductivity (SC), assuming that the sys- 
tem thermalizes. The obvious interest is that it is far 
easier to realize SC in attractive systems with s-wave 
pairing than in repulsive ones where the SC gap has to 
be anisotropic, so that the ac drive may provide an ad- 
vantageous as well as novel avenue for SC. However, a 
simple, sudden quench would heat the system to, e.g., 
\T cS /J oS \ = 1.83 ± 0.02 for A = 4,0 = 2tt, and U = 1 
(Fig. [4]) , which is much larger than the critical temper- 
ature for the attractive Hubbard model {T c /J ~ 0.1 for 
U ~ — 2 [HI]). S° the crucial question is whether we can 
avoid such a heating. 

A smooth change in U c s from positive to negative 
might seem desirable, but this is unfortunately impossi- 
ble, since \U e g/U\ has to cross a singularity at |j7o(»4)| = 
0. Instead, we propose here a "multistep quench" (an 
example is shown in the left panel of Fig. [5]). (i) We start 
from a weak interaction, (ii) This is followed by a sudden 
quench (to avoid the above problem) to an attractive but 
still weak interaction. The reason we start from a small 
U is that a larger jump of J7 e ff in a sudden quench tends 
to cause a higher temperature, (iii) After the sudden 
quench we employ a smooth ramp to amplify U c g . 

In the right panel of Fig. Owe plot T c s(t) for one exam- 
ple of the multistep U quench (U e g = 0.2 — > —0.5 — » —2), 
which is estimated by equating the total energy of the 
system at each rescaled t with the one for an equilibrium 



system with temperature T . For the example displayed 
here |T eff (t = 5)/J off | = 0.080 ± 0.004 with an initial 
T/J = 0.05, which accomplishes a temperature lower 
than T c /J ~ 0.1. 

Finally we mention the stability and experimental fea- 
sibility of the population inversion. Once the negative- 
T distribution is realized, it does not relax to a lower 
energy state due to energy conservation as long as the 
system is isolated from the environment with no energy 
dissipation. If fermions are coupled to other degrees of 
freedom, the distribution should start to collapse to the 
normal one with a positive T with a decay rate deter- 
mined by the strength of the coupling to the environ- 
ment. One good candidate that can avoid this difficulty 
is a system of cold fcrmionic atoms trapped in an opti- 
cal lattice, which is a virtually ideal, isolated single-band 
system [5j|. There, a conversion from attractive to re- 
pulsive interactions, which is equally feasible, may also 
be interesting. For electron systems, possible candidates 
are superlattices or arrays of nanostructures designed to 
realize a single-band well separated from other bands. 
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by the Japan Society for the Promotion of Science. 
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